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ABSTRACT 

Reactive melt infiltration is viewed as a promising means 
of achieving near-net shape manufacturing with quick 
processing time and at low cost. Since the reactants and 
products are, in general, of varying density, overall 
conservation of mass dictates that there is a force related to 
chemical conversion which can directly influence infiltration 
behavior. In effect, the driving pressure forces may compete 
with the forces from chemical conversion, affecting the 
advancement of the front. We have developed a two- 
dimensional numerical code to examine these effects, using 
reaction-formed silicon carbide as a model system for this 
process. We have examined a range of initial porosities, pore 
radii, and reaction rates in order to investigate their effects on 
infiltration dynamics. 

NOMENCLATURE 

Roman letters 

k r reaction rate constant 

MW ( species molecular weight 

p pressure 

r pore radius 

R ( rate of production/consumption of species t 
t ei time at which infiltration is complete 
r time at which chemical conversion is complete 

; 3 time at which infiltration enters third, post-reaction- 

dominated regime 

u D macroscopic Darcy (front) velocity 
v f species molar specific volume 

V volume 

>y location of the infiltration front 

Greek letters 

£ intrinsic infiltration rate 


K 

permeability 


molar consumption rate 

<t>C 

species volume fraction, V, jV 

<J> 

porosity 


volume fraction of solid, <X> 5 = 1 - 4> 

o 

surface tension 

0 

wetting angle 

V 

fluid viscosity 

S u bsc r i ots/su oe rsc ri Dts 

c 

capillary 

f 

front 

i 

Index denoting species 


BACKGROUND 

In reactive melt infiltration (RMI), a porous medium, or 
"preform”, is first created, which is easy to form into a desired 
shape. A liquid infiltrant is introduced into the preform through 
a pressure gradient caused by capillary action or by applied 
pressure (see, e.g., Kochendoerfer, 1991; Hillig, 1994a, b; 
Haggerty and Chiang, 1990). The chemical conversion 
between the infiltrant and the preform results in a final product, 
which may be fully dense and is of the same dimensions as 
the original preform (Singh and Behrendt, 1992). 

In comparison to other methods of composite 
manufacture, RMI is quick and efficient. Chemical vapor 
infiftration requires days or weeks to complete the chemical 
conversion process, but RMI needs only minutes or hours. 
Post-process machining is limited, making this process 
attractive for the manufacture of brittle and expensive 
ceramics. 

Reaction-formed silicon carbide is chosen as a model 
system for this study. The production of monolithic silicon 
carbide system is appealing because the chemistry and 
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preform geometry are reasonably well characterized. We will 
focus on infiltration that is driven solely by capillary pressure. 
The relevant thermophysical properties of the chemical 
constituents are given in Table 1 . 

Table 1. THermophysical properties of silicon, carbon, 
and silicon carbide. 


Property 

Si 

c 

SiC 

P (kg'm 3 ) 

2400 

1500 

3120 

H (kg/m s) 

5X10 -4 



MW (kg/kmol) 

28.1 

12.0 

40.1 

7 (m 3 /kmol) 

0.0117 

0.008 

0.0128 


Although the absolute viscosity of liquid silicon, p, is a 
weak function of temperature (Turovskii and Ivanova, 1974), 
we will neglect any temperature dependence. Note that the 
specific molar volume, v ( , varies substantially among the 
species. In fact, the volume occupied by a mole of SiC is 
approximately 35% smaller than that of a mole each of Si and 
C. Based on overall mass conversion, Singh and Behrendt 
computed that the critical porosity, <P C , for this process is 
0.358. For initial porosities below <D C , the flow will choke 
prematurely because there simply is not enough volume 
available to complete the conversion. If O 0 > O c and all of the 
carbon is converted in a fixed preform volume, the resulting 
ceramic will be comprised of SiC and Si. In this study, we will 
investigate O 0 in the range of 0,40-0.60, well above the 
critical porosity. 

In order to model the preform, we have drawn on the 
extensive work of Behrendt and Singh on the melt infiltration of 
silicon and silicon alloys into microporous carbon (see 
Behrendt, 1986; Behrendt and Singh, 1993, 1994; Singh and 
Behrendt, 1992, 1994a, b, 1995). Their work was also the 
basis for other experimental and analytical studies 
(Sangsuwan et al., 1998; Datta et al., 1996; Messner and 
Chiang, 1990). The amorphous, glassy carbon preforms were 
uniform and isotropic with good connectivity (Singh and 
Behrendt, 1992, 1994b; Sangsuwan et al., 1998). Moreover, 
they could be consistently reproduced by Behrendt and 
Singh’s process with pore diameters ranging from 0.1 -3.0 pm 
and initial porosities from 0.33-0.64. We will examine pore 
radii ranging from 0.10-0.20 pm and initial porosities from 
0.40-0.60. 

Infiltration takes place in an evacuated furnace. Without 
the application of external pressure, the infiltration process is 
driven by capillary action. Consequently, the infiltration rate is 
dependent upon the wetting characteristics of the infiltrant 
upon the preform surface, which is in turn sensitive to local 
temperature and surface and melt impurities. At 1703 K, the 
contact angle, 0, of liquid silicon on carbon is reported to be 
38° by Li and Hausner (1991). (See also Whalen and 
Anderson, 1974; Turovskii and Ivanova, 1974.) We will use a 
value of 0.725 N/m for the surface tension, a (Brandes, 1983), 
although other values have been reported (Minnear, 1982). 


We will assume that the chemistry is reaction-controlled 
rather than diffusion-dominated, based on the work of Messner 
and Chiang (1990) and Fitzer and Gadow (1984) so that the 
pore radii close by: 

r(f) = r 0 -k r t } . (1) 

where r 0 is the initial pore radius, k r is the reaction rate, and 
r y is the local reaction time, i.e., the amount of time that has 
elapsed since the reaction began in this region. Since the 
temperature variation within the preform for this set of 
parameters is on the order of 100°C (Nelson, 1998; Einset, 
1993), we will use a constant value for k r . Messner and 
Chiang estimated the reaction rate to be 4x1 0" 8 m/s, but 
bounded the estimate from 1 0" 8 to 1 0” 7 m/s. Alloyed melts may 
decrease reaction rates as well as eliminate the residual 
silicon phase (Singh and Behrendt, 1992, 1994a,b, 1995). 
Accordingly, we will investigate reaction rates in the range of 
10” 9 to 10- 7 m/s. In reality, the reaction rate is not a tunable 
parameter, and the effect on infiltration dynamics will vary if 
the chemistry is modified. 

GOVERNING EQUATIONS 

A schematic of the two-dimensional numerical model of 
reactive melt infiltration is sketched below in Fig. 1. Liquid Si is 
introduced at the bottom of the preform and wicks up into the 
preform due to capillary action. 



The general form for the conservation of mass is given for 
each of the species by: 

= t = Si, C, SiC (2) 

ot 
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( 7 ) 


where R t is the rate of production/destruction of species 
and <p ( is the species volume fraction, i.e., the fraction of the 
volume occupied by species L The species volume fraction 
must obey the constraints 0<<^<1 and < 1. In any 
region that is completely filled, X t 0< =1 * Since overall mass 
must be conserved, the conservation equation for the total 
mass in the system is subject to the constraint X/ Rt = 0 ■ The 
liquid phase, comprised of a single component, i.e., Si, is the 
only phase with a nonzero velocity. 

We will model this process as a zeroth-order chemical 
reaction. For stoichiometric conversion, this exothermic 
reaction is given by Si + C-> SiC , so that the local production 
of species can be expressed as: 


R t = firf MW f 


(3) 


3*5 

dt 


* 




r o 




J 


For details, see Nelson, 1998; Nelson and Colella, 1999a,b. 
Recall that r(t) = r 0 -k r t r , where r 0 is the initial pore radius 
and t r is the relevant reaction time. Combining Eqs. (5) and 
(7) gives the molar consumption rate as: 

n= - 1- — (8) 

v SiC v c r o V J o ) 

At the pore level, the convective time scale is much longer 
than the diffusive time scale. Therefore, we can use d’Arcy’s 
law to describe the motion of the front: 


where = +1 for the product and -1 for the reactants, r| is the 
number of moles consumed per unit time, and MW t is the 
species molecular weight. We note that none of the species 
can be updated to a negative value. Since we are provided 
with an abundant source of Si, the amount of C available 
locally will limit the chemicaf conversion. Therefore, when all of 
the C has been consumed locally, the chemical reaction 
comes to a halt in that region. 

We will make the critical assumption that the overall 
volume of the preform remains constant (Singh and Behrendt, 
1992). Using the above assumptions, the overall conservation 
of mass can be written as: 

^ * (05/ 1 ~ 5 / ) = r f ( y sic ~ ^si “ T c ) (4) 

where v, =MW, jp ( is the specific molar volume. 

The expression for rt is found by relating the geometric 
description of the pore closure to that obtained from the 
conservation of mass. The volume fraction of the solid phase 
is <X> 5 = <p c + <j> SiC . From mass conservation, the change in solid 
volume fraction with respect to time is: 

d<P s _ d<t> c d(psic 

dt dt dt (5) 

= *7 ( v 5rC “ V c) 

If we envision the preform porosity to be comprised of a 
bundle of parallel capillary tubes, then the permeability for this 
model is found by: 



(6) 


where r is the pore radius and O is the porosity ( 0> 5 = 1 - <X> ). 
Experimental data (Datta et al. p 1996) indicates that this is an 
adequate model for this preform. Comparing the cross- 
sectional area of pore space initially and at time n, we find 
that: 


11 D = 



(9) 


where u D is the front velocity, k is the permeability of the 
porous medium, \x is the absolute viscosity of the fluid, and p 
is the pressure. For this system, k is a function of space and 
time. Furthermore, any body force term is negligible in 
comparison to the pressure force provided by capillary action. 

We can relate the fluid velocity, ii Si , to the front velocity, 
U D , by the Dupuit-Forcheimer relation, i.e., u D = <t> Si u Si . We do 
not explicitly know u D or the pressure field; however, our 
analysis of continuity provided a relationship for the 
divergence of velocity in Eq. (4). Combining the divergence of 
d’Arcy’s law, Eq. (9), with Eq. (4), we arrive at a steady 
Poisson problem for pressure that stems from the 
conservation of momentum and of mass: 


r K ^ 

“Vp 
^ J 


- *7 { V SiC ' 


( 10 ) 


Equation (10) identifies the mechanism by which capillary 
forces and chemical conversion can enter into competition in 
the infiltration process. For this system, the right-hand side of 
Eq. (10) is at minimum zero. When it is nonzero, the pressure 
profiles within the preform exhibit some curvature, resulting in 
a nonuniform velocity field. 

For this system, the ambient pressure is negligible. We 
apply a homogeneous Dirichlet boundary condition at the 
lower boundary; at the infiltration front, the pressure is equal to 
the negative of the capillary pressure, -p c . Using the Young- 
Laplace equation, we can estimate p c at the front as a function 
of surface tension, c, pore radius, r , and the contact angle, 0: 


Pc = 


2cr cos 0 
r 


( 11 ) 


Depending on the initial pore radius, p c varied from about 
5.94x1 0 6 to 1 .88x1 0 7 Pa. As in the experimental setup of 
Sangsuwan et al. (1998), we will assume that the preform is 
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coated with, e.g., boron nitride, to inhibit mass transport across 
the sidewalls. This ensures that the liquid silicon from the bath 
will not wick up the sidewalls, and that no fluid will escape from 
the system through the sidewalls. Therefore, homogeneous 
Neumann boundary conditions are applied to the pressure at 
the side boundaries. 

After the pressure solve, we can directly calculate the 
velocity from the pressure field by using d'Arcy’s law, Eq. (9). 
The normal component of velocity at the sidewalls is zero, due 
to the homogenous Neumann boundary condition used on the 
pressure field. Although we would expect some tangential 
velocity slip in the immediate vicinity of the sidewalls, these 
effects are confined to distances on the order of a pore 
diameter (Kaviany, 1991). Therefore, we will assume that, at 
the left and right boundaries, dv/dx = 0 . In this set of 
numerical simulations, the infiltration process is essentially 
one-dimensional, as in experimental work of a similar process 
(Einset, 1993). This can be expected due to uniform initial 
porosity and pore size, the choice of boundary conditions on 
pressure and velocity, and the simplicity of the chemistry and 
permeability models. 

The derivation of the energy equation, as weii as the full 
details of the above derivations, can be found in Nelson (1998) 
or Nelson and Colella (2000b). 

NUMERICAL PROCEDURE 

The equations derived above are used to develop a finite- 
volume algorithm for reactive melt infiltration. The discrete 
formulation in two dimensions Is discussed in Nelson (1998) 
and Nelson and Colella (2000a). Briefly, the macroscopic 
velocity is used to advance the front on a fixed mesh using a 
front-tracking technique developed by Chem and Colella 
(1988), coupled with a volume-of-fluid method and front 
reconstruction. Chemical reaction is allowed to occur based on 
the local species volume fractions and local reaction time. The 
updated species volume fractions are used to form a new 
permeability field. Next, a Poisson problem is solved on the 
infiltrated region only to find the new pressure field, using the 
technique of Johansen and Colella (1998). Finally, d’Arcy’s law 
is used to find the new velocity field. The process is repeated 
until all of the C in the preform is consumed or the final 
simulation time is reached. 

The mesh size was 32x128 on a 0.01x0.04 m domain so 
that the mesh spacing Ax = Ay = 0.0003125 m. A description 
of the discrete equations, as well as the effects of grid 
resolution and other sources of numerical error are discussed 
in Nelson and Colella (2000a) and Nelson (1998). 

RESULTS 

In this section, we will explore the effects of varying 
reaction rate, initial porosity and initial pore radius. 


Effect of reaction rate 

With the initial porosity fixed at 0.50 and the initial pore 
radius at 0.17 pm, the infiltration curves for reaction rates of 0, 
10~ 9 , 10" 8 , and 10~ 7 m/s are shown below in Fig. 2. The 
capillary pressure corresponding to this case is 6.99x1 0 6 Pa. 
Fig. 2(a) presents the location of the front at the centerline, 
>y(0 , while Fig. 2(b) shows >y 2 (r). 

In the absence of chemical reaction, y f {t) exhibits a 
square-root dependence, and the slope of >y 2 (r), called the 
intrinsic infiltration rate (Hillig, 1994b) and denoted as is 
constant. The front reaches the upper boundary of the 
preform, located at 0.04 m, at =15.51 s. 



Figure 2. Centerline infiltration profiles for O> 0 =0.50; 
r 0 =0.17 |im; k r =0 (1), 10‘ 9 (2), 10' 8 (3), and 10' 7 (4) m/s. 

(a) >y(0; (b) y) 2 (if 

For the chemical system of Si/C/SiC , the infiltration 
process is retarded by the chemical reaction. For the low 
reaction rate, k r = 10” 9 m/s, the infiltration curves closely match 
the nonreacting case for the first couple of seconds. After that 
point, the front begins to lag behind the nonreacting case. 
Fig. 2(b) shows that that \ is not constant over this time 
interval. In contrast, y f 2 for A r =10" 7 m/s appears quite linear 
after half a second. For k r - 10 -8 m/s, there is a short period of 
under one second during which the curves closely follow the 
nonreacting profiles. Up to about 8 seconds, the slope of y f 2 
is continuously changing. After that time, £ appears constant, 
nearly matching that of k r =1 0“ 7 m/s. 

What explains this behavior? The infiltration moves 
through a sequence of flow regimes in which the forces 
induced by capillary action and by chemical reaction are in 
competition (Nelson, 1998; Nelson and Colella, 2000a, b). 
Initially, the pressure forces dominate as the infiltrant 
impulsively moves into the preform, indicated through the 
close correspondence of yj with the nonreacting case. This 
can be seen in approximately the first two seconds of £ r =10" 9 
m/s and the first second of k r =1 0 _B m/s. 

In the second regime, the forcing term provided by the 
chemical reaction flares up and wanes. Within the preform, 
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Figure 3. Centerline profiles of (a) pressure; (b) velocity; (c) permeability; and (d) <p c for <t> 0 =0.50, ^ =0.17 pm, and 
k r =10' 8 m/s. Profiles are presented in approximately one-second increments for t =1-5 seconds, and in ten-second 

increments thereafter. 


the centerline pressure profiles become progressively more 
curved (see Eq. (13)), leading to decreased flow velocities 
(see Eq. (9)), as shown in Figs. 3(a) and (b) for ^=10"® m/s. 
Behind the front, the permeability, k, Fig. 3(c), and <t > c , 
Fig. 3(d), steadily decrease. Previous work (Nelson, 1998, 
Nelson and Colella, 2000a, b) indicates that the centerline 
thermal maximum lags well behind the infiltration front in this 
regime. Furthermore, the maximum operating temperatures 
were observed in this regime, particularly with insulated 
boundaries nearby. 

When all of the available C has been consumed near the 
entry to the preform, Fig. 3(d) at /=7.99 s, the pressure 
relaxes to a more linear profile except in the immediate vicinity 
of the front, Fig. 3(a). Correspondingly, the velocity profile 
becomes constant over part of the infiltrated region, Fig. 3(b). 
Consequently, {t) exhibits the characteristic bump in the 
infiltration profiles of Fig. 2, e.g., in k r = 10" 8 m/s between 4 
and 8 seconds. This behavior can be observed in experimental 
work by Einset (1996). 

For ^=10“® m/s, the infiltration is still early in the 
reaction-dominated regime at 15 seconds, while for£ r =10“ 7 
m/s, the infiltration has passed through the reaction-dominated 


regime by about one-half second. As will be seen throughout 
this work, the temporal and physical extent of the reaction- 
dominated regime decreases with increasing reaction rate. 

After 10 seconds, the infiltration process for fc r =10”* m/s 
enters a third regime, in which the reaction-generated forces 
are confined to a small region behind the front and the 
infiltration dynamics are largely set by the permeability in the 
fully reacted zone and the capillary pressure. From 
10 seconds to complete infiltration at 63 seconds, y f {t) 
resembles a square-root profile, and \ can be approximated by 
a constant. The reaction zone becomes progressively more 
localized to the vicinity of the front. For example, the reacting 
zone at 63 seconds is about 0.00296 m, appearing over eight 
computational cells. As a result, the thermal maximum is very 
close to the infiltration front (Nelson, 1998). 

Although the curve for the highest reaction rate in 
Fig. 2(b) appears linear after one second, notice that, unlike 
the nonreacting case, the curve does not intersect the origin. 
Examination of the initial time interval more closely indicates 
that the infiltration simply proceeds through the reaction- 
dominated regime more rapidly than the other cases, resulting 
in the offset. If we were to extend the temporal window for 
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Figure 4. Effect of <J> 0 on infiltration profiles for r 0 = 0.17 pm and 4> 0 =0.40 (1), 0.45 (2), 0.50 (3), 0.55 (4) and 0.60 (5): 
k r =0 (black), 10 9 (blue), 10^ (green), and 10 7 m/s (red): (a) y f vs. t\Q>)y f 2 (t) vs. /. 


k r = 10” 9 m/s to later times (and with a longer preform length, in 
most cases), we would also see the characteristic bump 
signaling the ascendance and conclusion of the reaction- 
dominated regime, but occurring at a retarded pace. 

When fully reacted, the final value for <p Sl is 0.197 and for 
(j) sic * it is 0.803. 

Effect of initial porosity 

To examine the effect of varying porosity, infiltration 
profiles are shown in Fig. 4 for £ r =0, 10~ 9 , 10 -8 , and 10~ 7 m/s 
with <& 0 =0.40, 0.45, 0.50, 0.55 and 0.60, holding r 0 fixed at 
0.17 pm. The capillary pressure corresponding to this r 0 is 
6.99x10 6 Pa. In this time interval, the infiltration profiles for 
k r =0 and 10” 9 m/s collapse into two curves. The increased 
permeability caused by increasing the porosity by 50% does 
not significantly impact the infiltration in the pressure- 
dominated regime and the early stages of the reaction- 
dominated regime. 

At increased reaction rate, we can still see an early period 
during which the curves for the varying porosity are 
indistinguishable. However, in the subsequent reaction- 
dominated regime, Fig. 4 indicates that the reaction-generated 
forces can greatly slow down the infiltration, as seen most 
dramatically for k r - m/s at the lowest porosities. In fact, 
the reaction-generated forces can theoretically stop the flow 
for a time. The fanning effect seen for k r = ICf 8 and 10“ 7 m/s 
can be explained as follows: Since carbon is the only species 
present in the preform initially, <& Q is related to the amount of 
C present, i.e., <f>c 0 = 1 - <J>o- This in effect places an upper 
bound on the amount of chemical conversion that is possible 
over a given volume. Once the carbon in the initial reaction- 
dominated region is consumed, the dynamics of the infiltration 
proceed at a different pace. At constant k r , examination of 


the centerline profiles reveal that the fanning occurs soon after 
the pressure profile relaxes back from its most distorted 
shape. These trends indicate that the reaction-dominated 
regime becomes larger in spatial and temporal extent with 
decreasing porosity and reaction rate. For the fully reacted 
zone, the final </> SlC was 0.96, 0.88, 0.80, 0.72, and 0.64 for 
# 0 =0.40, 0.45, 0.50, 0.55 and 0.60, respectively. 

The times required for complete infiltration and complete 
conversion, t ei and t e r , are presented in Fig. 5. We did not 
carry out the infiltration process to complete conversion for 
k r =10~ 7 m/s due to numerical instability; this is being 



Figure 5. Influence of <J> 0 on required time for complete 
infiltration (unfilled symbols) and conversion (filled 
symbols) for r 0 = 0.17 pm and A r =0, III;10 d , O;and10 a 
m/s, ♦. 
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addressed in a modified algorithm (Nelson and Colella, 
2000a). 

Notice that the most efficient reaction rate for infiltrating 
the preform, i.e., the lowest k f . , does not imply the shortest 
conversion time, particularly at low initial porosities. Also note 
that, as the reaction rate is increased at fixed porosity, the 
difference between t ec and r eJ decreases, it is zero for 
<£ 0 =0.40, £,.=10“ 8 m/s, indicating that the carbon behind the 
front is essentially consumed as fast as the front can move. 



(a) 





Figure 6. Influence of <X> 0 in the post-reaction-dominated 
regime: (a) time of entry into this regime; and (b) intrinsic 
infiltration rate. 


The time at which the infiltration moved from the reaction- 
dominated regime to the third post-reaction-dominated regime, 
r 3 , shown in Fig. 6(a), was estimated from infiltration profiles. 
The intrinsic infiltration rate, in this regime is presented in 
Fig. 6(b). 


For Jt r =0, ^I.OxlO" 4 m/s over the range of porosity. For 
k } . =10" 8 m/s and 10 -7 m/s, we can see that q increases 
linearly with porosity. We might expect this, because the 
permeability is proportional to the porosity, while capillary 
pressure remains fixed. When chemical reaction occurs, the 
entry into the post-reaction-dominated regime is retarded as 
porosity decreases, reflecting the increased influence of the 
reaction-generated forces on the infiltration dynamics. 
However, once in the post-reaction-dominated regime, the 
permeability in the fully reacted zone behind the front is more 
influential on the infiltration dynamics than the reaction- 
generated forces. 


Effect of initial pore radius 

Next, we will examine the influence of the initial pore 
radius on the infiltration process, holding the porosity fixed at 
0 o =O.5O and varying the reaction rat e,fc r , to 0, ICf 9 , 1CT 8 , 
10“ 7 m/s and the pore radius to 0.10, 0.14, 0.17 and 0.20 pm. 
The capillary pressure corresponding to these pore radii 
ranges from 5.94x1 0 6 Pa to 1.88x10 7 Pa (see Eq. (11)). For 
clarity, the infiltration curves are presented in sets by reaction 
rate in Fig. 7 for the first 30 seconds of the infiltration. 

Decreasing pore radius retards the infiltration, as 
observed by Einset (1996) experimentally. Ignoring the 
chemical reaction for the moment, this would be expected 
since the capillary pressure is inversely proportional to r 0 , 
while the permeability is proportional to r 2 . Unlike the 
parameterization in porosity, the curves for k r =0 and 10~ 9 m/s 
do not collapse onto one another in Figs. 7(a) and (b). We also 
see the effects of chemistry in the curves for k r =1 0 -8 and 1 Cf 7 
m/s in Figs. 7(c) and (d). 

The curves for fc,.=1CT 9 m/s (Fig. 7(b)) are nearly 
equivalent to the corresponding nonreacting case for the first 
second or two. All of the preforms at this reaction rate are fully 
saturated by 30 seconds except for /q= 0.10 pm. The intrinsic 
infiltration rate, shown in Fig. 7(f) is not constant for any of 
the pore radii over this time interval. 

The characteristic bump in the infiltration profile due to the 
presence of the reaction-dominated regime is seen most 
clearly in Fig. 7(c) for k r =10“ s m/s under 10 seconds. The 
corresponding slope of y f 2 (t) in Fig. 7(g) is not linear during 
this time. The infiltration has moved beyond the reaction- 
dominated regime by 10 seconds in Fig. 7(c). This occurs 
more quickly for the higher reaction rate in less than a second 
in Fig. 7(d). From Figs. 7(g) and (h), we can see that £ 
appears constant once the infiltration process has left the 
reaction-dominated regime. 
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time (s) time (s) time (s) time (s) 


(a) (b) (c) (d) 



time (s) time (s) time (s) time (s) 


(e) (f) (g) (h) 

Figure 7. Effect of ^ on infiltration profiles for O 0 =0.50; r 0 =0.10 (1), 0.14 (2), 0.17 (3) and 0.20 (4) pm: (a) >y (r), =0.0; (b) 

y f (t\ k r =10' 9 m/s; (c) y f (t\ k r =10' 8 m/s; and (d) y f (t\ k r =10' 7 m/s; (e) k r - 0.0; (f) y f 2 {t\ A r =10' 9 m/s 5 (9) 

y f 2 {t\ k r = 10* m/s; and (h) y f 2 (t), k r = 10 7 m/s. 


The dependence of the total infiltration and conversion 
times on pore radius is shown in Fig. 8. A wide gap is once 
again apparent between t eJ and t e r for £ r =1Cf 9 m/s (unfilled 
and filled circles). For the lowest pore radius and reaction rate 
(filled circle, left of arrow), the infiltration process had moved 
well beyond the reaction-dominated regime when the upper 
boundary was encountered. On the other hand, much of the 
conversion took place in a fully saturated preform for the larger 
pore radii (filled circles, right of arrow). Recall that, for a fully 
saturated preform, the pore radius decreases linearly with 
time, Eq. (1), to complete conversion. The trend of increasing 
t e r with increasing radius at high pore radii therefore indicates 
that, under these conditions, the absence of the front reduces 
the effectiveness of the conversion. 

Figure 9 shows the influence of pore radius on the start of 
the post-reaction-dominated regime and on the intrinsic 
infiltration rate. Fig. 9(a) indicates that decreasing pore radius 
reduces the temporal extent of the reaction-dominated regime, 
an effect which is most pronounced at the lower reaction rate. 
Unlike the porosity-dependent parameterization, % is not 
constant for k,.=0 (squares), but rather increases with 
increasing pore radius, due to increasing permeability. As 
before, we should expect this because the capillary pressure is 
inversely proportional to r 0 , while the permeability is 
proportional to r 0 2 . Fig. 9(b) indicates that 4 is linearly 


proportional to r 0 and is. not affected greatly by the reaction 
rate. These results imply that the infiltration in the post- 
reaction-dominated regime is governed primarily by the 
permeability of the preform in the reacted zone behind the 
front and by the capillary pressure. 



r 0 (x10' 7 m) 


Figure 8. Influence of r Q on required time for complete 
infiltration (unfilled symbols) and complete conversion 
(filled symbols) for k, =0, □; 10 -9 , •; and 10" 8 m/s, ♦. 
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Figure 9. Influence of r 0 in the post-reaction-dominated 
regime: (a) time of entry into this regime; and (b) intrinsic 
infiltration rate. 

CONCLUSIONS 

We have examined the infiltration of liquid Si into a 
constant-volume microporous C preform to form a monolithic 
SiC ceramic composite for the range of O 0 =0.40-0.60, 
r 0 =0.10-0.20 pm, and k r - 0, 10~ 9 , 10 , 10“ 7 m/s. The initial 
pore radius and porosity were uniform over the entire preform. 

We have seen that the infiltration dynamics are dependent 
on reaction rate, initial porosity and initial pore radius. 
Significantly, it is also a function of the flow regime arising from 
the competition between pressure and reaction forces. We did 
not find a simple correlation for the time and distance 
associated with the reaction-dominated regime. However, 
numerical simulation can be used to predict these details. 

Initially, the pressure forces dominated as the infiltrant 
impulsively entered the preform. Later, the presence of the 
chemical reaction retarded the infiltration due to a competition 
between the pressure forces and the volumetric forces due to 
the chemical conversion. When the reactants in the entry to 
the preform were consumed, the reaction zone was confined 
to a small region behind the front. Furthermore, the infiltration 


rate increased, although it was significantly smaller than that of 
the corresponding nonreacting case. 

In the first pressure-dominated regime and the early 
stages of the second reaction-dominated regime, a porosity 
increase of 50% had a negligible effect on the infiltration 
profiles, as seen in the data for k,.= 0 and 10~ 9 m/s. 

The reaction-generated forces may impede or enhance 
the flow, depending on the density variation among the 
products and reactants. For the Si/C/SiC system with a dense 
product, the conversion was shown to retard the flow. The 
temporal and physical extent of the reaction-dominated regime 
decreased with increasing reaction rate, decreasing pore 
radius and increasing porosity. In this regime, the velocity 
within the preform was retarded significantly due to the 
curvature of the pressure profiles, causing a characteristic 
bump in the infiltration profiles. Theoretically, the flow can 
become stationary at the front; this is to be distinguished from 
choking in that this is a temporary phenomenon. For k r = ICf 8 
m/s, the reaction-dominated regime was in evidence from 
about 2 to 10 seconds, while it was completed by one second 
for k r = 10" 7 m/s. The intrinsic infiltration rate was not constant 
in this regime, indicating that the location of the front did not 
have a square-root dependence on time. Aside from the 
implications on infiltration dynamics, the maximum process 
temperature was found, in previous work, to occur in this 
regime. The reaction-dominated regime ended soon after the 
carbon in the initial entry zone was completely consumed. 

Beyond the reaction-dominated regime, the infiltration 
settled into a simpler pattern. The reaction zone was localized 
to an area just behind the front. The intrinsic infiltration rate 
was constant and proportional to the initial porosity and to the 
initial pore radius. Although the permeability is proportional to 
the square of the pore radius, the capillary pressure is 
inversely proportional to the pore radius, accounting for this 
behavior. The intrinsic infiltration rate was not greatly 
influenced by the magnitude of the reaction rate in this regime. 
We therefore deduce that the infiltration was primarily affected 
by the capillary pressure at the front and the permeability in 
the fully reacted zone behind the front. Consequently, the 
intrinsic infiltration rate can be characterized by a constant 
beyond the reaction-dominated regime, although the line 
representing £ t must be offset from the origin to coincide with 
the data. 

For this system, numerical simulation of the parameter 
space provides the best prediction of optimal parameters for 
minimum conversion time. For an actual process, we should 
also include preform length as a parameter, as that will affect 
whether or not the process leaves the reaction-dominated 
regime. For k r = 10" 9 m/s, the infiltration was completed during 
the early stages of the reaction-dominated regime, except at 
the smallest pore radii. If the entire preform reached saturation 
during this time, the overall conversion time was seen to 
increase with increasing pore radius, since the evolving pore 
radius decreased strictly linearly with time, as set by the 
chemistry of the system. In contrast, if the infiltration moved 
into the post-reaction-dominated regime, as seen at lower pore 
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radius, the conversion time decreased with increasing pore 
radius. The presence of the advancing front, with its 
accompanying complex infiltration dynamics, can either 
improve or reduce the efficiency of the conversion process, 
depending on the initial parameters. If the infiltration was 
allowed to proceed through the sequence of infiltration 
regimes, the overall conversion time decreased with 
increasing porosity, pore radius and reaction rate. 

This sort of analysis could provide insight into the design 
of functionally graded materials. For example, local 
modification of the pore radius or porosity in the post-reaction- 
dominated regime could be used to speed up the infiltration 
rate. Such a modification in the reaction-dominated regime 
could have major implications on the efficiency of the process. 
The advantage of a numerical simulation is that it could range 
over the desired parameter space with minimal effort, allowing 
for more precise optimization. 

References 

Behrendt, D.R. “Porous silicon carbide as a matrix for 
ceramic composites.” NASA TM 88837 (1986). 

Behrendt, D.R. and M. Singh. “Effect of carbon preform 
pore volume and infiltrants on the composition of reaction- 
formed silicon carbide materials.” J Mat Synthesis and 
Processing 2.1 17-123 (1994). 

Behrendt, D.R. and M. Singh. “Theoretical considerations 
for reaction-formed silicon carbide.” NASA TM 106414 (1993). 

Brandes, E.A., ed. Smithefls Metals Reference Book, 6 th 
ed. Butterworths (1983). 

Chern, l-L and P. Colella. “A conservative front tracking 
method for hyperbolic conservation laws”. Unpublished paper 
(1988). 

Datta, S.K., N. Simhai, S.N. Tewari, J.E. Gatica, and M. 
Singh. “Permeability of microporous carbon preforms." Met 
Trans A 27 A: 3669-74 (1996). 

Einset, E. O. Private communication (1993). 

Einset, E. O. “Capillary infiltration rates into porous media 
with applications to Silcomp processing.” J Am Ceram Soc 
79:333-338(1996). 

Fitzer, E. and R. Gadow. "Investigations of the reactivity of 
different carbons with liquid silicon." In Proceedings of 
International symposium on ceramic components for engines. 
Japan, 1983. KTK Scientific Publishers, Tokyo, Japan (1984). 

Haggerty, J.S. and Y.M. Chiang. “Reaction-based 
processing methods for ceramics and composites.” Ceram 
Eng and Sci Proc 11: 757-781 (1990). 

Hillig, W.B. “Making ceramic composites by melt 
infiltration.” Am Ceram Soc Bull 73 56-62 (1994). 

Hillig, W.B. “Melt infiltration approach to ceramic matrix 
composites.” Am Ceram Soc, Comm 71:C-96 to C-99 (1994). 

Hillig, W.B., R.L. Mehan, C.R. Morelock, VJ. DeCarlo and 
W. Laskow. “Silicon/silicon carbide composites.” Amer Ceram 
Soc Bull 54: 1054-56 (1975). 


Johansen, H. and Colella, P. “A Cartesian grid embedded 
boundary method for Poisson’s equation on irregular 
domains.” J Comp Phys 147: 60-85 (1998). 

Kaviany, M. “ Principles of heat transfer in porous media . ” 
Springer-Verlag (1991). 

Kochendoerfer, R. “Liquid silicon infiltration: a fast and 
low-cost CMC manufacturing process." In Composites: 
Proceedings of the 8th International Conference on Composite 
Materials. Held in Honolulu, HI, July 15-19, 1991, pp. 23-F-1 to 
23-F-9.3 (1991). 

Li, J.G. and H. Hausner. ‘Wettability of silicon carbide by 
gold, germanium and silicon.” J Mat Sci Ltrs 10. 1275-76 
(1991). 

Messner, R. and Y.M. Chiang. "Liquid-phase reaction- 
bonding of silicon carbide using alloyed Si-Mo melts." J Amer 
Ceram Soc 73. 1 1 93 ( 1 990) . 

Minnear, W.P. "Interfacial energies in the Si/SiC system 
and the Si+C reaction." J Amer Ceram Soc 65: C10-C11 
(1982). 

Nelson, E.S. “A numerical algorithm for the simulation of 
reactive melt infiltration.” Ph.D. Thesis. University of California 
at Berkeley (1998). 

Nelson, E.S. and P. Colella. “A numerical algorithm for 
reactive melt infiltration.” To be submitted to J Comp Phys 
(2000a). 

Nelson, E.S. and P. Colella. “Infiltration dynamics in 
reactive melt infiltration.” To be submitted to J Amer Ceram 
Soc (2000b). 

Nelson, E.S. and P. Colella. Unpublished work (2000c). 

Sangsuwan, P., S.N. Tewari, J.E. Gatica, M. Singh and R. 
Dickerson. “Reactive infiltration of silicon melt through 
microporous amorphous carbon preforms.” Met Trans 
30B$ 33-944 (1998). 

Singh, M. and D.R. Behrendt. “Reactive melt infiltration of 
silicon-molybdenum alloys into microporous carbon preforms.” 
Mat Sci Eng A 194: 193-200 (1995). 

Singh. M. and D.R. Behrendt. “Microstructure and 
mechanical properties of reaction-formed silicon carbide 
(RSFC) ceramics.” Mat Sci Eng A 187: 183-1 87 (1994). 

Singh, M. and D.R. Behrendt. “Reactive melt infiltration of 
silicon-niobium alloys in microporous carbons.” J Mat Res 
7:1701-1708(1994). 

Singh. M. and D.R. Behrendt. "Studies on the reactive 
melt infiltration of silicon and silicon-molybdenum alloys in 
porous carbon." NASA TM 105860 (1992). 

Turovskii, B.M. and U. Ivanova. "Temperature 
dependence of the viscosity of fused silicon." Inorganic 
Materials 70:2108-2111 (1974). 

Whalen, T.J. and AT. Anderson. “Wetting of SiC, ShNa, 
and carbon by Si and binary Si alloys." J Am Ceram Soc 
58: 396-399(1974). 


NASA/TM— 2000-209802 


10 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-01 B8). Washlnglon, DC 20503. 


AGENCY USE ONLY {Leave blank) 


TITLE AND SUBTITLE 


REPORT DATE 

March 2000 


3. REPORTTYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Parametric Study of Reactive Melt Infiltration 


6. AUTHOR(S) 


Emily S. Nelson and Phillip Colella 


WU-101-43-0C-00 


PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-12107 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 2000-209802 


11. SUPPLEMENTARY NOTES 


Prepared for the 1999 International Mechanical Engineering Congress and Exposition sponsored by the American 
Society of Mechanical Engineers, Nashville, Tennessee, November 14-19, 1999. Emily S. Nelson, NASA Glenn Re- 
search Center, and Phillip Colella, Lawrence Berkeley National Laboratory, M-S SOD- 133, Berkeley, California 94720. 
Responsible person, Emily S. Nelson, organization code 6727, (216) 433-3268. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Categories: 24, 27 and 34 


12b. DISTRIBUTION CODE 


Distribution: Nonstandard 


This publication is available from the NASA Center for Aerospace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

Reactive melt infiltration is viewed as a promising means of achieving near-net shape manufacturing with quick process- 
ing time and at low cost. Since the reactants and products are, in general, of varying density, overall conservation of mass 
dictates that there is a force related to chemical conversion which can directly influence infiltration behavior. In effect, the 
driving pressure forces may compete with the forces from chemical conversion, affecting the advancement of the front. 

We have developed a two-dimensional numerical code to examine these effects, using reaction-formed silicon carbide as a 
model system for this process. We have examined a range of initial porosities, pore radii, and reaction rates in order to 
investigate their effects on infiltration dynamics. 


14. SUBJECT TERMS 

Infiltration; Porosity; Composite materials; Silicon carbides; Porous materials; Reacting 
flow; Fluid dynamics; Manufacturing; Production engineering; Permeability 


15. NUMBER OF PAGES 


17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

20. LIMITATION OF ABSTRACT 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 


Unclassified 

Unclassified 

Unclassified 



Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 























